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We present the exact solution to the linearized Maxwell equations in space-time slightly curved by 
a gravitational wave. We show that in general, even dealing with a first-order theory in the strength 
of the gravitational field, the solution can not be written as the sum of the flat space-time one and a 
weak perturbation due to the external field. Such an impossibility arises when either the frequency 
of the gravitational wave is too low or too high with respect to the one of the electromagnetic 
field. We also provide an application of the solution to the case of an electromagnetic field bounced 
between two parallel conducting planes. 
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I. INTRODUCTION 



O . 

The problem on the propagation of electromagnetic fields in a curved space-time has been widely studied all through 
the years. Einstein himself calculated the deflection of a light ray by the gravitational field of the Sun. This problem 
can be simply approached by assuming light rays as trajectories of a zero-mass particle, for which ds 2 = 0. Moreover 
the same approach has been used to calculate the response of an interferometer to an impinging gravitational wave [Q . 
, A different approach to the problem is considering directly the Maxwell equations in curved space-time. Two kinds 
of phenomena can be described by this method: propagation of electromagnetic radiation in a non-radiative gravita- 
tional background [0-|| (see also Ref. § for an exhaustive treatise on the subject) and in the field of a gravitational 
wave. Not only are they interesting problems themselves but also have important experimental applications (see for 
instance ||). 

As for the propagation in a gravitational wave field, some authors have solved the problem by using the formalism of 
qq ' geometrical optics; therefore these results are valid under the assumption that the gravitational wave is slowly varying 
with respect to the electromagnetic field ||-[nj. Others have solved the equations by splitting the electromagnetic 
tensor (or the four-vector potential) in a sum of two terms: the unperturbed solution and a correction term [p2|-p0[ . 
However, as it has been shown in Ref. [ jn] , the splitting approach turns the problem of the propagation of a free 
electromagnetic field in a gravitational background into a different one, that is the generation of a field by a given 
far\' current. On this ground many doubts could be raised as to the validity of this procedure (a detailed discussion may 
be found in Ref. Q). 

Therefore we see that a satisfactory solution that only assumes the smallness of the gravitational wave amplitude 
regardless of its shape or frequency is lacking. In this paper we provide exact solution to the problem by solving 
Maxwell equations for the 4-potcntial (De Rham equations) in the framework of the linearized general relativity (Sec. 
II). In other words the solution we have achieved is correct as long as the linearized De Rham equations ( |2.3| ) describe 
properly the physical problem. We express our result as an integral which is the curved space-time version of the 
usual Fourier integral. In this way any solution can be obtained by a proper choice of the Fourier-like coefficients. 
Afterwards we calculate the expressions of the electromagnetic tensor components (Sec. III). Finally Sec. IV is 
devoted to the application of our solution to a simple case in order to compare it to the one obtained through the 
"splitting" procedure. 

Besides we find that the main effect of a low-frequency gravitational wave on an electromagnetic field bounced 
between two parallel conducting walls is ohmic power loss. 



II. EXACT SOLUTION TO THE DE MAXWELL EQUATIONS 



The homogeneous de Rham equations describe the propagation of electromagnetic field in an external curved 
space-time. Under the Lorentz condition, they could be written as: 
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A*? a = 
A°' = 



(2.1) 



(throughout this paper the notation and conventions are as in plf exept for indices p, r, s which take values 1 and 
2). Let us assume that the curvature of space-time is due to a weak plane gravitational wave, propagating in the x 3 
direction of a given transverse traceless (TT) reference frame plf . In this case the metric tensor is: 



\hnu\ < 1, 



(2.2) 



where r)^ v is the Minkowski metric tensor with positive signature and h^{x 3 — x°) is the perturbation to the flat 
space-time metric. In the chosen reference frame the only non -vanishing components are /in = — h 22 = h+ and 
hi2 = h 2 i = h x . Neglecting second order terms in h^, Eqs. (2.1) become {h^ v = rj^ a h au ): 



Aw 



A\ = 



h^A^p + (h^J + h^ a - hj>>») A a fj = 



(2.3) 



Since the first equation of the above system for fi = is the same as \i = 3, we can choose a particular Lorentz gauge 
in wh ich A 3 = A . This greatly sim plifies the problem, because this way the equations are partially decoupled [see 
Eqs. (|2.5|) -(2.8) later on]. Eqs. (2.3) are a set of differential equations with variable coefficients, which depend only 



This fact leads us to perform the following coordinate transformation: 



(x 3 



x°)/V2 



(x 3 +x°)/V2 



(2.4) 



With this new set of variables, system (2.3) becomes (<9 M means 



A° = Q; 



cM 1 



DA P + d h p s d 3 A s = 
DA 3 - d h r s d s A r = 
DA = 



d 2 A 2 



d 3 A 3 = 0. 



(2.5) 
(2.6) 
(2.7) 
(2.8) 



where 



d = (did! + d 2 d 2 + 2d 3 d - h rs d r d s ) , 

p,r,s = 1,2, and A^ — A v are the new components of the 4- vector potential. We observe that in this new 
reference frame the electromagnetic gauge condition is the so-called Coulomb gauge O]. Hereinafter only the real 
part of equations is to be retained. First of all we solve Eq. (2.5); through Fourier transform in (X — > (X 1 , X 2 , X 3 ), 
A — > (Ai, A2, A3)) 



A P (X,X") = 



0\ _ / e »A fc X 



e 



a p (\,X°)d 3 \, 



(2.9) 



one gets: 



21X3 Jx° ~ ^ + x ^ + XrXshr s (^ X ° 



dh p 

a p + iX 3 ^a r =0 



(2.10) 



Setting: 



i p = 1 (X°; \)B p (V2X°; A) 



(2.11) 



with 



7 (X U ; A) = cxp 



(Af + X 2 )X° - X r X s h^ rs (V2X°) 



2*Aa 



(2.12) 
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(h^ 1 ) rs stands for any primitive of h rs ) Eq. (2.1C) gives: 



dB p 1 dh p 



lB q = 0. 



dX° 2 dX° 

The solution to the above equation could be expressed as a function series: 



(2.13) 



oo (-iy 



K r Jn;w) &p(A) 



* r P (n;w) = ftUw p ( n -l;w) dw 



(2.14) 



pv ' 7 J dw 



where b p (\) depends only on A. As one can easily see, H r p (n;w) are of same order as h n , where h is the order of 
magnitude of the gravitational wave amplitude. A sufficient condition for the series to converge is that H r p (n; w) < 
1; this is always satisfied under the assumption of lineariz ed th eory, in which h -C 1. Therefore, through this 
approximation, we can keep only the first term in the series ( 2.14 ). In this case 



B r (72 A ; A) = b r {\) - h r p (y/2 A ) b p (X). 



As for the A 3 component, the Coulomb gauge [Eq. (2.8)] gives: 

o Ai a 1 + A2 a 2 



A, 



(2.15) 



(2.16) 



which automatically fulfills Eq. (2.6). 

Having solved the system of equations [2.1)— ( pjq ) in the A M reference frame [Eqs. (2.4)], we proceed to write the 
solution in the x a coordinates. In the case of a flat space-time, i.e. = 0, kj wou d be the Fourier variables 

conjugated to coordinates. In this way we have that A 3 = (k 3 ± k)/\/2 (k = \Jk\ + k 2 + fc|). One gets (k — > 
(ki,k 2> k 3 )): 



A p = I B p (x 3 -x°;k)e l 



A 3 = A" = / k s B s (x 3 -x°;k)e 



kx 
k x 



/3+(fc)e l(fe4+,/ ' +) + (3-{k)e-' l{kt -^ ) 



d 3 k 



k + k 3 



,i (kt+i> 



+ ) P (fc) -Hkt-i/,-) 



k - fes 



d 3 k 



where 



± _ k r k s h^ rs {x 3 -x°) 
^ ~ 2 (fc ± k 3 ) 



(2.17) 
(2.18) 

(2.19) 



and P ± (k) are generic functions whose actual form is determined by assigning suitable initial or boundary conditions. 
The solution given in Eqs. (2.17) and ( 2.18| ) onl y ass umes the smallness of and does not require any further 
hypothesis about its rate of variation. From Eq. (2.19) it follows that: 



x 



where \ is the order of magnitude for a typical frequency of gravitational wave. As no assumption was made over x, 
ip is not necessarily a first order term in h. Therefore, in the general case, it is not possible to expand the exponential 
terms in the solution to the lowest order in -0, opposite to the case of the function series (2.14), in which there is no 
term proportional to 1/x- 

If the gravitational wave is linearly polarized then: 



h+(x 3 -x°) = H+f(x 3 -x°), 
where H + and H x are constants. We can set: 

H + = H cos 2a, 



h x (x 3 -x°) = H x f(x 3 -x°) 



H sin 2a. 



(2.20) 



(2.21) 
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In this case the solution to Eq. ( |2.13 ) great ly simplifies. Actually it is possible to express it as a function rather than 
an infinite series of functions [see Eq. ( 2. 14 )]. Using standard methods the solution can be written: 



B p = b p f(x ' 
where b p and A are the eigenvectors and eigenvalues of the following system: 

H p q b q = Xb p . 

The eigenvalues are A = H and A = —H while the corresponding eigenvectors are: 

b + = (cos a, sin a) and b = (sin a, — cos a) 
respectively. Therefore the solutions read 



(2.22) 
(2.23) 
(2.24) 



Bl=b p + e- ! t f M and B p _ = b p _e+^ f ^ . 



(2.25) 



III. THE ELECTROMAGNETIC TENSOR 

Having found the solution to the de Rham equations, it is straightforward to determine the components of the 
electromagnetic tensor by means of 



5 1 — A — A 



We explicitly have: 



(3.1) 



F u 



,kx 



where (dot means derivation with respect to a; 3 — x°): 

± ..A. i 

J 03 



fop 



ft, 



fp3 



hp k s 



1 \ h h h 7 

±k [b p + ^h p r b r )Tb p 1 ^Pf- 1 , 

! 2 r J 2 (k ± fc 3 ) k ± fc 3 



1 



1 



kp k s 



1 \ h h h rs 

b p + -h p b r ] =F fr p r s I 

+ 2 r J T ° 2 (jfc ± fc 3 ) A; ± k 3 



b s - h s r b r 



b s - - h s r b r 



h p r b r 



-\h p r b r 



(3.2) 

(3.3) 
(3.4) 
(3.5) 
(3.6) 



M axw ell equations in vacuum, F 11 ^ = 0, in the linear approximation of the chosen TT gauge reference frame [see 
Eq. Q] read: 



( v pa tff* ~ ^ a h vf} - rf fi h pa ) F aP , v + hf" T) vfi F a/3 = 0, 



(3.7) 



wher e terms of t he o rder of h 2 and x^ 2 are negligible. Actually the electromagnetic field components [Eqs. (3.2)- 
(3.6)] satisfy Eq . (|3.7| ) with the same degree of precision. The homogeneous Maxwell equations are identically fulfilled 
because of Eq. (|3.1| ). 

Therefore our solution solves the problem of the propagation of an electromagnetic field in a gravitational wave 
background with the only assumption of a small perturbation to a flat metric, thus generalizing the result of fll|| to 
any gravitational wave frequency. 
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IV. APPLICATION 



As an example let us consider the case of an electromagnetic field bounced between two conducting parallel walls a 
distance L apart. In the absence of gravitational waves, the physical features of the problem we are concerned with do 
not depend on the specific orientation of the apparatus. This is due to the isotropy of flat space-time. However this 
isotropy breaks down in the presence of a radiative gravitational field propagating along a given direction. Because of 
this a nisot ropy, elect romagnetic waves propagating in opposite directions are in general affected in different ways [see 



Eqs. (2. IE) and (3.2)]. A consequence of this circumstance is the impossibility of obtaining stationary solutions. Only 
if the propagation of the electromagnetic field lies on the plane perpendicular to the propagation of the gravitational 
wave, a stationary solution does exist. In fact, only in th is case, the perturbation to the phase of the electromagnetic 



wave is the same for opposite directions [see Eqs. ( 2.19| ) and (3.2)]. The peculiarity of this plane also arises if we 



notice that two different points lying on it undergo the same gravitational perturbation. 

In Subsec. A we consider the last case and compare our result with the already existing solutions obtained by 
splitting of the electromagnetic field |l2|,^0| (see Sec. I). We find out that the latter approach lies in the framework 
of linearized theory only in the definite range of gravitational wave frequencies given by the relation: 

k 

fc/i<X<T (4.1) 

h 

This is different from what was pointed out in Ref. p^| , where such a condition was mandatorily needed for the 
linearized theory to be valid. Indeed this is merely the condition for the validity of the splitting procedure. We think 
that this circumstance is due to the expansion in power of the gravitational wave amplitude of the electromagnetic 
phase, expansion that is needed by the splitting method. This kind of procedure (expansion of a phase) is not in 
general correct, as it is shown, for instance, by the long experience in perturbative methods in both classical and 
relativistic celestial mechanics (e.g. calculation of the advance of the perihelion of Mercury plfl). 

In Subsec. B we show, in a particular case, how the non stationariness of the solution induces dissipative effects on 
the conducting planes. 



A. Comparison with the splitting procedure 



As a first case we assume the walls to be orthogonal to the x 1 -axis, placed at a; 1 = and x 1 — L, respectively. Let 
us consider the particular solution to the de Rham equations that, when h nv is switched off, describes a stationary 
electromagnetic field propagating along the x 1 direction with a linear polarization along the :r 2 -axis. For the sake of 
simplicity we also consider a linearly polarized gravitational wave for which H x — (see Sec. II). This solution is: 



Fm — 



E 



F 



23 



—2 B n sin (knX 1 
y^gg» sin{k n x 1 ) 



(f - Hf) sin k n x° 



HK 



2 k n 



fsin{knX o + EhL f (-D 



cos k„x 







Hk n 



K-i) 



1 12 



= J2 2S n ( 1 - ^ / ) COS (Kx 1 ) COS ( k n X° + ^/ 



Hk r , 



K-i) 



(4.2) 

(4.3) 
(4.4) 



where k n = and n is an integer. We notice that the solution is still stationary in the a^-direction. One can 
see that the gravitational wave induces a new component in the electromagnetic tensor (namely F23). For a direct 
comparison to Ref. [Q we set f{x 3 — x°) — — cosx(a; 3 ~ x°), and B n — B n 5 nn . Performing a formal expansion in H 7 
the electromagnetic field components are found to be a sum of the flat space-time solution and a perturbation, that 
is: 



F 



(0) F„„ + (%, 



//// 



where: 



(0) ^02 = -2 B ri sin {k n x l ) sin {k n x°) 
(o)F 23 = 

(°)Fi2 = 2 B fl cos {k n x 1 ) cos {k n x°) 



(4.5) 



(4.6) 
(4.7) 
(4.8) 
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are the "unperturbed" solutions, while 



«F 02 = -2HB n sin^x 1 ) 
W F 23 = -HB n sin (fcfix 1 

^F 12 =HB n COS^X 1 ) 



cosx(x 3 — x ) sm(k fi x 1 ) — 



1 / kn 



+ — I smx(x 3 ~ x°) cos(knX°) 



cosx(a; 3 — a; ) sin (k n x°) — — sinx(z — t) cos (k n x°) 



cosx(a; 3 — a; ) cos(/c fi a; 1 )H — — sinx(z — t) sin(fcfjX°) 



(4.9) 
(4.10) 
(4.11) 



would represent the perturbation due to the gravitational wave. Actually from the above exp ress ions one immediately 
achieves that the formal expansion of the exact solution is meaningless when condition (hi) does not hold true. 
In fact, in this case the "perturbation" terms are even greater than the "unperturbed" ones. On the other hand if 
Eq. ([hi]) does hold, the splitting approach Jl2|,^(| yields a result that coincides with the formal expansion in H of our 
solution, but in a rotated reference frame. 

It is important to notice that the solution expressed by Eqs. (4 1 2)-(h4) can be cast in the form 

A = Y / a n (x 7 x 3 -x°)A n (x) 7 (4.12) 

n 

where A stands for any electromagnetic tensor component, and A n are the same as in a flat space-time. The effect 
of the gravitational wave results in a modulation of the a n coefficients, and not in a creation of any new mode A n > 
(n' ^ n) even if \ — \k n ± k n ,\. 



B. Dissipative effects 



As a second case of application we consider the walls to be orthogonal to the bisector to the x —x plane, placed in 
x 1 = —x 3 and x 1 — \[2L — x 3 , respectively. In a similar way to what we have done in the previous case, we consider a 
solution of the de Rham equation describing, when h uv vanishes, a stationary electromagnetic field propagating along 
the x 1 = x 3 direction with a linear polarization along the x 2 -axis. We set h x — again. One has: 



n 

H 



x 1 + x 3 



Hf (-i) 



B n Hf 
V2 



/cos k n x + 



k n 



(1-Hf)sm(k n x° 



k n Hf(~ 



cos 



71 



x 1 + x 3 



Hf(~ 



cos I k n x° 



F 23 =^-\/2B n (l-fT/) cos 

n 

+ B n H sin 



71 



x 1 + x 3 - 



Hf 



k n 

.71 

2 

k n 

71 



X + x — 



Hf (-D 



X + X 



f sin ^k n x° - 
Hf (-D 



cos yk n x° 
2 

cos I k„x 



— cos I k n x v H 



(4.13) 



(4.14) 



o 



(4.15) 



As in the previous case, the gravitational wave does not create any new modes of oscillation A n i [see Eq. ( 4. 12 )]. 

The most interesting feature of the solution is that the parallel component to the surface of the electric field does 
not vanish on the walls any more. This results in a gravitationally induced ohmic power loss within the conductors 
which is to be added to usual losses (such as those due to thermal vibrations of the conducting planes) already existing 
in flat space-time and not depending on the orientation of the apparatus. Consequently if the former kind of loss is 
at least comparable to the latter one it would be possible, in principle, to infer the presence of a gravitational wave 
by measuring the time decay of a resonant cavity. 

We now proceed to estimate the gravitationally induced power loss per unit surface when 



Hk n 



> H. 



(4.16) 



G 



This assumption leads to interesting results. In order to study dissipative currents induced by the gravitational wave 
on the walls, we only need to retain the biggest part of the field; it is noticeable that all terms proportional to H are 
negligible when Eq. ( 4.16| ) is fulfilled. One gets {B n = 5 n nBn): 



-2B„ sin * 



F 23 = -V2B n cos * - 



2V2 



sin ykfix + 
cos ( knX° 



knHf 



(-1) 



2:! 



(4.17) 

(4.18) 
(4.19) 



where = if x 1 = —x 3 and ^ = kn L when x 1 = —x 3 + \/2L. If we set U — ^p-, where U is the electromagnetic 
energy and da stands for the surface element parallel to the walls, then the ohmic loss is given by p3| : 



du_ 

~dt 



,3 J? 2 

c E \\ 
8n S lu 



(4.20) 



where c is the speed of light (in the following we drop the convention c = 1, hence x° — ct and fl — \ c stands for 
the gravitational wave frequency), Eu is the electric field component parallel to the surface, 5 is the skin depth, and 
lu = ckn the electromagnetic frequency; factor 2 takes into account the fact that there are two walls. Therefore the 
equation describing the rate of energy loss per unit surface is: 



1 dU 
U ~dt 



2 c 2 
LSu 



sin 



knHf(~V 

2V2 



(4.21) 



where an average over — has been performed. The solution to the above equation can be easily obtained assuming 
X x 3 ~0: 



U{t) = Uq exp 



2 c 2 
L5 uj 



sin 



2V2 



dt' 



(4.22) 



If w e set f(x) = cos (%x) and make the further assumption (whose validity must be checked at the end) Ot < 1, 
Eq. ([4.22Q becomes: 



The time after which the energy falls by a factor 1/e (time decay) is: 



fl 2 t 3 



(4.23) 



3 LSl 



2 n 2 c 2 



2V2~n 



(4.24) 



If the gravitational wave source considered is a binary pulsar as, for instance the PSR 1913 + 16 pit], then we find: 



t = 5 (5L)% 10 4 sec - 50 sec < ^ - 5 10 3 sec 



(4.25) 



where we have set uj = 10 15 radsec 1 and {5 L) 



10- 



7 2 /3 



(e. 



Unfortunately this time lies many order 



of magnitude above the usual — non gravitationally induced — optical resonant cavity decay times. 
However if we consider a coalescing compact binary [ ^6| we get: 



sin(2fir) _ 8SLn 3 



c 2 ojH 2 



(4.26) 



where Eq. ( 4.22| ) has been used. Analysis of this equation, by taking mi — TO2 — 3 m Q , shows that the shortest decay 
time is obtained when the time to coalescing ranges within 2 10 6 sec ~ 23 d. In this case r — 0.8 sec, to be compared 
with the decay times of the order 1-10 msec in a cavity. Therefore we see that this time is long, as compared to 
presently-measured decay times, but not so long to rule out the possibility to measure it in the future. 
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V. CONCLUSION 



We have found the exact solution to the linearized Maxwell equations for a gravitational background described by 
a plane gravitational wave. We have not made any assumption over the shape or frequency of the gravitational wave. 
The only approximation used is the linearized Einstein theory, within which Eqs. (3.7) hold. 

As an application we have evaluated the electromagnetic field bounced between tw o co nducting walls. A comparison 
with known approximated results has been performed. We have found that, if Eq. (4.1) is met, the formal expansion 
in H of our solution reduces to the app roxim ated one. Besides we have found out that the gravitational wave does 
not cause any new mode A n i [see Eq. ( 4.12 )] to appear in the cavity. Furthermore we have evaluated a possible 
measurable effect of the interaction with low frequency gravitational waves. We have shown that in general an ohmic 
loss on the walls is produced; this results in a decrease of the stored electromagnetic energy. 
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